function [error,dependency]=calc_error(chi,covar)

	indx=find(covar);
	[n,m]=size(covar);
	redcov1=covar(indx);
	nn=sqrt(length(redcov1));
	redcov=reshape(redcov1,nn,nn);
	iredcov=inv(redcov);
	codia=diag(redcov);
	icodia=diag(iredcov);
	error=sqrt(chi*diag(covar));
	dependency=1-1./(codia.*icodia);
